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Abstract 

The role of the rigidity of a peptide chain in its equilibrium dynamics is in- 
vestigated within a realistic model with stringent microscopically derived coupling 
interaction potential and effective on-site potential. The coupling interaction char- 
acterizing the chain rigidity and the effective on-site potentials are calculated for 
three main types of protein secondary structure. The coupling interaction is found to 
be surprisingly weak for all of them but different in character: repulsive for a-helix 
and anti-parallel /3-sheet structures and attractive for parallel /3-sheet structure. 
The effective on-site potential is found to be a hard one for a-helix and anti-parallel 
/3-sheet and a soft one for parallel /3-sheet. In all three types of protein secondary 
structures a stable zig-zag shape discrete breather (DB) associated with the os- 
cillations of torsional (dihedral) angles can exist due to weakness of the coupling 
interaction. However, since the absorption of far infrared radiation (IR) by proteins 
is known to require the existence of rather long chains of hydrogen bonds that takes 
place only in a-helicies, then one can conclude that the excitation of a DB in such 
a way is possible in a-helix and seems to be hardly possible in /3-sheet structures. 
The interpretation of the recent experiments of Xie et al. on far IR laser pulse 
spectroscopy of proteins is suggested. The frequency of a DB in an the a-helix is 
obtained in the region of 115 cm -1 in accordance with these experiments. 



1 Introduction 



Nonlinearity opens new perspectives for solving the problem of localization and 
storage of energy by biomolecules [1], [2], [3]. This problem is especially urgent 
within the context of enzymatic catalysis. How does an enzyme store and 
utilize the energy released at substrate binding? The answer to this question 
is very poorly understood not only at quantitative but even at qualitative level. 
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By now it is clear that dynamical contributions from linear modes into the 
mechanism of enzyme action are negligible because of high friction in such a 
condensed matter as protein interior leading to fast dissipation of energy [4]. In 
contrast nonlinear systems provide us with examples of surprising behavior: 
the higher the friction is the slower the dissipation of energy proceeds [5]! 
Though the generosity of this striking result is under question one can hope 
that since proteins are very rich structures for exhibiting nonlinear phenomena 
then some of them may play a functional role in enzymatic reactions. Among 
nonlinear phenomena in biomolecules the so called discrete breather (DB) or 
else intrinsic localized mode was put in the forefront during the last decade 
and even overshadowed such a fashionable concept as soliton. However most 
results on DBs were obtained within unrealistic models with toy potentials. 
That is why the progress in application of this new paradigm to proteins 
requires taking into account detailed information on protein dynamics if one 
wish to pass from toy models to realistic ones. 

Last years are marked by noticeable progress in revealing subtleties of protein 
dynamics gained by infra-red (IR) spectroscopy [6], [7], [8], [9], [10], [11], nor- 
mal mode analysis [12], [13], [14] and molecular-dynamics simulations (see the 
recent reviews [15], [16] and refs. therein). On the one hand time-resolved 2D 
measurements data suggest that in a-helicies there are distributions of the tor- 
sional (dihedral) angles (ip and ip) around their average values with a width of 
20 degrees and that the conformation fluctuates on a time scale of picoseconds 
[7] confirming previous molecular-dynamics studies. On the other hand the ex- 
periments of Xie et al. [10], [11] demonstrate that very long-lived oscillations 
can exist in the a-helix proteins (myoglobin, bacteriorhodopsin) but not in the 
/3-sheet one (photoactive yellow protein) when excited by a far IR laser pulse. 
The oscillations are observed in the Amide-I band (u ~ 1500 cm" 1 ) [10] and 
more vividly in the low-frequency band u ~ 100 cm" 1 [11]. The absorption of 
radiation by proteins in the latter case is commonly attributed to the excita- 
tion of longitudinal and transverse modes of one-dimensional hydrogen bonded 
chains in a-helices [17]. The normal mode analysis of myoglobin [12], [13], [14] 
testifies that the anharmonic decay rate of higher frequency localized normal 
modes, calculated by perturbation theory, is typically nearly independent of 
temperature, consistent with results of pump-probe studies on myoglobin. In 
fact, the long-lived vibrations in the Amide-I region of myoglobin [10] are not 
all that long-lived, about 15 ps, so that this is really a predictable result. 
Indeed, Leitner [14] and Yu and Leitner [12] have presented calculations of 
anharmonic decay rates in the Amide-I region of myoglobin, revealing that 
vibrational lifetimes of order 1-10 ps are quite likely and unsurprising. On the 
other hand, the results of [11] remain a puzzle. In this study, long-lived exci- 
tations of over 500 ps are found near 115 cm" 1 . Though there are no direct 
comparisons for bacteriorhodopsin, the anharmonic decay rates in myoglobin 
near 100 cm" 1 computed by Yu and Leitner are two orders of magnitude faster 
than this. 
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According to the conjecture of the authors of [10], [11] their result suggests 
that a breather can be created in an a-helix. Another interpretation of these 
experiment within the framework of Davydov's soliton concept is also sug- 
gested [18]. The latter utilizes the possibility of longitudinal mode excitation 
at the absorption of far IR radiation by proteins. However this interpretation 
encounters with difficulties at the point emphasized in [10] that a-helices in 
myoglobin or bacteriorhodopsin are too short for such soliton could efficiently 
propagate in them. In this chapter the hypothesis is put forward that the 
long-lived oscillations in a-helicies are the fluctuations of the torsional angles 
of peptide groups around their equilibrium values. The aim of this chapter is 
to derive the possibility of the existence of a DB in polypeptide chain from its 
microscopic model and to interpret the phenomena observed the papers [10], 
[11]. Our approach emphasizes another possibility as opposed to the paper 
[18]: the excitation of transverse mode of one-dimensional hydrogen bonded 
chains in a-helices at the absorption of far IR radiation by proteins is actually 
that of oscillations of the peptide groups around their equilibrium positions 
in peptide chain (the oxygen and nitrogen atoms oscillate transverse to the 
chain of peptide groups). This type of motion is shown to enable the exis- 
tence of DBs in a-helices. Thus the DB in the peptide chain of a-helices can 
provide sustaining the long-lived oscillations excited by the absorption in the 
low-frequency band cu ~ 100 cmT 1 . 

DBs or else intrinsic localized modes (time periodic spatially localized oscilla- 
tions with significant amplitudes of several units in a chain of weakly coupled 
non-linear oscillators while others are at rest or oscillate with negligible ampli- 
tudes) discovered by Sievers and Takeno in 1988 [23] and proved by MacKay 
and Aubry [24] to be structurally stable are well understood and commonly 
appreciated as generic phenomenon in nature at present (see [25] and refs. 
therein). The key ingredient for the existence of a DB is the requirement 
of weak coupling interaction between adjacent nonlinear oscillators [24], [25]. 
DBs have become a new and very fruitful paradigm in nonlinear physics. The 
concept of DB has been applying to DNA dynamics for a long time [19], [1], 
[2], [3], [20], [21]. Its application to protein dynamics is damped by a point of 
view (see e.g. [3]) that proteins are less capable than DNA to exhibiting DBs 
because they are much less regular structures. However there are fragments of 
secondary structures in proteins that are highly regular. In [26] a model for a- 
helix protein is considered that exhibits DB existence in the chains of peptide 
groups connected by hydrogen bonds in the spirit of Davydov's model for a 
soliton. In this chapter a model is constructed that pursues the same goal but 
for the peptide groups connected by C - C bonds in secondary structures of 
the backbone. Another distinction from [26] and from all other paper on DBs 
in biomolecules is that we deal with realistic effective on-site potentials and 
coupling interaction potentials rather than with model toy ones. These poten- 
tials are derived from a stringent microscopic model of the polypeptide chain 
and thus the present work seems to be a step to ab initio calculation of DBs in 
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biomolecules. There are few analytical models treating non-linear mechanics 
of peptide chain backbone [2], [27], [28], [29]. The reason for transparent lack 
of activity on analytical modeling of peptide chain backbone mechanics seems 
to stem from meager direct experimental evidence for nonlinear excitations in 
it. Nevertheless the experiment [11] can be considered as that and motivates 
the present attempt to construct a simple and tractable but at the same time 
microscopically stringent model of polypeptide backbone dynamics. 



2 The Hamiltonian of the model 

A schematic picture of a polypeptide chain with all designations used further 
is presented in Fig.l. The mutual orientation of two adjacent peptide groups 
is characterized by the torsional (dihedral) angles (pi and ip^. The angle ipi 
characterizes the rotation round the bond Cf — iVj and that ipi characterizes 
the rotation round the bond Cf — C i . Further we consider the equilibrium 
dynamics in protein secondary structures in which the equilibrium values of the 
torsional angles are the same for all peptide groups (they are listed at the end 
of the next Sec). We consider small deviations tpi(t) and ipi(t) of the torsional 
angles from their equilibrium values (ifi = (p® + (fii(t) and ipi = ip^ + ipiit)) 
with \ipi(t)\ < 20° and \ipi(t)\ < 20°. At such amplitudes of the deviations the 
hydrogen bonds confining the peptide group in the secondary structure are 
not broken [4]. Finally we consider a peculiar type of motion ipi_i(t) = <fi(t). 
The latter means that the peptide group rotates as a whole respective some 
effective axis a passing through the center of the C— N bond parallel to the 
bonds Cf - ^ and N i+1 - Cf +1 that are assumed to be approximately parallel 



(ZCf -Cl- N i+1 = 113° and ZC,' - N i+1 - Cf +1 = 123° [30]). This type 



of motion is stipulated by the fact that the peptide group is a planar rigid 
structure [4], [30]. For the sake of uniformity of designations we further denote 
Xi = 2ipi(t) = 2ipi_i(i) (see Fig. 2) and thus 



The moment of inertia of the peptide group relative to the axis a can be easily 
calculated and is / ~ 7.34* 10~ 39 g* cm 2 . The Hamiltonian of the polypeptide 
chain in our model with nearest neighbor interactions is 



Here Ui oc (xi) includes interactions defining the local potential of the peptide 
group (namely hydrogen bonds, the so-called torsional potentials and the van 



(1) 
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der Waals interaction of covalently non-bonded atoms of the peptide group 
with the atoms of adjacent side chains) defining its separate motion while 
U(xf,Xi + i) includes the coupling interactions (namely the van der Waals in- 
teraction of covalently non-bonded atoms of adjacent peptide groups i and i+1 
and their electrostatic interaction) intermixing the motions of these groups and 
leading to the rigidity of the peptide chain. It should be stressed that the lat- 
ter potential also contributes into the separate motion of the peptide groups. 
We can define the effective on-site potential for such motion as 

V eff (xi) = U loc ( Xi ) + U( Xi] 0) + U(0; Xi ) (3) 

and the coupling interaction potential as 

U(x t ; x i+1 ) = U vdw ( Xi ; x i+1 ) + U el (x t ; x i+1 ) (4) 



Thus the equation of motion is 

j d 2 Xi _ dUi oc ( X j) _ dU (xj-i\Xj) _ dU(xj; x i+1 ) 
dt 2 dxj dxj dxj 



In the following two sections we consider in details the functions Ui oc (xi) and 
U(xi ; Xi + i). At doing it we will freely pass back and forth between the variables 
according to the rule (1) which can be rewritten as 

A<p i = x i /2] A^i = x i+1 /2 (6) 



3 Coupling interaction defining the rigidity of a polypeptide chain 



Both interactions in (4) intermixing the motions of the adjacent peptide groups 
and contributing to the coupling interaction are described by the central po- 
tentials W^™(R mn (<pi\ ipi)) and W^ n (R mn ((pi; ipi)) between the atom A n of the 
i-th peptide group and the atom A m of the i+l-th one with R mn being the 
distance between the atoms A n and A m . The electrostatic potential is 

WL (Rmn) = ^ (7) 



where q m and q n are partial charges on the atoms A m and A n respectively 
(q(N) = -0.28e; q(H) = 0.28e; q(0) = -0.39e; q(C) = 0.39e where e = 
4.8* 10 -10 CGS [30]) and e is the dielectric constant which for protein interior 
should be better conceived as some adjustable parameter (ss 2 4- 10 with 
3.5 being commonly accepted value close to high frequency permeability of 
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peptides) [30], [4]. In some studies e is supposed to be solvent dependent and 
chosen, e.g., 4 for CCl 4 , 6-=-7 for CHCl 3 and 10 for H 2 [31], [32]. For the 
van der Waals potential one can choose any of the numerous forms suggested 
in the literature, e.g., the Lennard- Jones one (6-12) or the Buchinghem one 
(6-exp). In this chapter we use for numerical estimates the former one 

A R 

rirvdw/n \ ^mn ^>mn / fi s 

mn mn 



with the set of the well known parameters of Scott and Scheraga (other sets 
were also verified and found to give similar results). 



The distance R mn as a function of the angles fi and is 

R mn (tPi,A) = {P m + r n+ 

2p m ^n[sin #(cos7 m sin a n cosy?j — sin7 m coso; ri cosipi) — 

sin7 m sin ^(cos^j cosy?j cos^ — sin-^j sin tpi) — cos7 m cos a n cos^]} 1 ^ 2 (9) 



The interaction potential for both van der Waals and electrostatic interactions 
has the form 

( vdw 1 f vdw 1 

U\ el I(x i; x i+1 )=J2w\ el i(Rmn(x i ;x i+1 )) (10) 



where the summation in n is over atoms in the i-th peptide group and that in 
m is over atoms in the i+l-th one. 

In what follows we work with full realistic coupling interaction potential de- 
scribed above. However for making use of suggestive analogies with the known 
literature results obtained on toy models we calculate the so called coupling 
constant which is a key characteristic of the truncated coupling interaction 
potential. Expanding the potential we obtain that to the leading order in the 
terms Xi 1 and Xi+i <^ 1 the rigidity of the peptide chain is determined by 
the term of the potential which intermixes the motion of the adjacent peptide 
groups 

U mix (x i ;x i+1 )K(-K)x i x i+1 (11) 



where the coupling constant is 
—K = -(K vdw + K el ) = 



6 



d 2 U- dw ( Xl ;x l+1 ) + d 2 U d (Xj,x i+1 ) 

9Xidx i+ i \ Xi =0;x i+1 =0 OXidXi + i \ Xi =0;x i+1 =0 



The results of calculations of the contributions into the coupling constant for 
different types of protein secondary structures are as follows: 

a) , a-helix (right) y?° = -57°; ^° = -47° [30]:-K vdw « 4.24 • 10~ 15 erg; 
-K el ~ 2.4 • 10~ 15 erg at e = 3.5; —K « 6.24 • KT 15 erg at e = 3.5. 

b) . anti-parallel /5-sheet $ = -139°; V° = 135° [30]:-K vdw « -5.0 • 1(T 16 
erg; -fT e ' w 2.64 • 1(T 15 erg at e = 3.5; —K « 2.14 • 10~ 15 erg at e = 3.5. 

c) . parallel /3-sheet y?° = -119°; = 113° [30]: -K vdw « -8.3 • 10~ 15 erg; 
-fT e ' w 2.96 • 10~ 15 erg at e = 3.5; -tf w -5.34 • 10" 15 erg at e = 3.5. 



4 Local potential of the peptide group 

The local potential Ui oc (xi) is composed as follows 

Uioc(xi) = U sc (xi) + U hb (xi) + U tors (xi) (13) 

Here U sc (xi) is the energy of van der Waals interactions of the atoms of the 
i-th peptide group with those of the side chains Ri and Ri-i and also with the 
atoms H l a and H l ~ l , Uhb(xi) = U$ (xi)+uj^ (xj) is the energy of two hydrogen 
bonds of the i-th peptide group and U tors (xi) = Ut ors {^i) + f/^^j-i) is the 
energy of the torsional potentials for the rotation of the angles ipi and V'i-i- 
For the latter we take the usual form 

U& ra (<p i ) = E v (l + cos3<p i ) (14) 

f/ls(^) = ^(l + cos3^) (15) 

with E v «1 kkal/mol and «1 kkal/mol. In all further simulations we 
choose the side chain to be Alanine {Ri is Cp(H) 3 ) and initially find the value 

of the angle xi ( see Fig.l) to minimize the energy. For the uj^(xi) (j=l,2) we 
take the potential [31], [32] 

uj&\ Xi ) = D(l - exp[-n(r(xi) - r )]f - D (16) 

where n=3 A -1 , r =1.8 A, Dm 4 6 kkal/mol at e — 3.5 is an adjustable pa- 
rameter of the hydrogen bond energy with D=5 kkal/mol being a conventional 



7 



value [4] (in [32] D is assumed to be a decreasing value at increasing e with 
D=0.5 kkal/mol at e — 10) and r(x.j) is the current length of the hydrogen 
bond. The dependence of the latter on Xi requires some simple and straight- 
forward but very cumbersome trigonometry which we omit here to save room. 
It should be only mentioned that this dependence can be expressed via the 
sets of the adjustable parameters P, L, 77 and R, L, k (see Fig.l). Here L is 
the length of the hydrogen bond at equilibrium (xj=0) (in the literature the 
values from 1.8 Atill 2.5 Aare figured), 77 and k are the angles characterizing 
the extent of linearity of the hydrogen bonds and finally R and P are the 
distances from crossing of the perpendiculars of the Oj_i and Hi atoms on the 
axis a to their hydrogen bond partners respectively at equilibrium (xj=0). In 
the present model we assume the partners to have fixed positions. Since in the 
a-helix the partners come from third peptide groups along the chain from a 
given one this assumption imposes a restriction on the size of a DB which can 
be considered, namely the assumption can be justified for DBs comprising no 
more than 5 sites (peptide groups). The elimination of the restriction for the 
a-helix would actually mean taking into account the long range interaction 
Uhb(xi] Xi± 3 ) along the peptide chain. Such complication is left for future work 
and in this chapter we exemplify the existence of a DB in the a-helix by a 3 
- site one. For /3-sheets the restriction can be discarded because the partners 
usually belong to very distant peptide groups from a given one. 

The adjustable parameters are chosen to satisfy the following set of require- 
ments: 1. For all three types of protein secondary structure the effective on-site 
potential V e ff(xi) must have a minimum at equilibrium (xi = 0) because the 
structures are steady stable ones. 2. For the anti-parallel /9-sheet the angles 77 
and k are known to be very close to 180° [30]. 3. For the a-helix the spectro- 
scopic frequency 



must be 115 cm -1 for the model to give interpretation of the experiment of [11]. 
Here A is a wavelength, c is the light speed and the dash denotes a derivative 
in Xi. These requirements are satisfied, e.g., at the following values of the 
common adjustable parameters: e = 3.5, D=5.7 kkal/mol, E v = 1 kkal/mol 
and — 1 kkal/mol and the following values of the particular adjustable 
parameters: 



a) , a-helix (right) y?° = -57°; $ = -47° [30]: L=2.026 A, R=3.68 A, P=3.34 
A, 77 = 170°, k = 170°. 

b) . anti-parallel /3-sheet y?° = -139°; $ = 135° [30]: L=2 A, R=3.61 A, 
P=2.82 A, 77 = 180°, k = 180°. 





(17) 



sp 



2lTC 
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c). parallel /5-sheet y?° = -119°; $ = 113° [30]: L=1.887 A, R=3.55 A, P=3.03 
A, r] = 170°, k = 170°. 



It seems plausible that the found values of the parameters are not unique and 
some variation of them compatible with the requirements mentioned above is 
possible but we could find no reasonable ways to it. However full exploration 
of the parameter space of the present model is a daunting task and we can 
not exclude the possibility that the found values of the parameters differ from 
reality Nevertheless considering them as most reasonable and coinciding with 
the conventional figures we further restrict ourselves by dealing with them 
only 



5 Results 

In what follows we use dimensionless variables. We define the frequency 



and measure time in the units of u; 1 so that the dimensionless time is 



Also we measure energy in the units of luJ 2 and denote the dimensionless 
coupling constant p 



The results of the investigation of the effective on-site potential and the cou- 
pling interaction potential can be summarized as follows: 



a) , a-helix ((1/A) sp =114.75 cm" 1 ) : p = -0.00193; V^fo = 0) > (for 
x < 0.3 the expansion is « 0.5a; 2 + 0.017a; 3 -0.235a; 4 + 0.6a; 5 -1.936a; 6 - 
3.8a; 7 + 8.437x 8 + 14.455a; 9 + 0.094a; 10 - 42.98a; 11 + 47.6a; 12 + ...). The effective 
on-site potential is presented in Fig.3. 

b) . anti-parallel /3-sheet ((1/A) sp =84.34 cm" 1 ) : p = -0.00115; = °) > 
(for x < 0.3 the expansion is « 0.5a; 2 + 0.004a; 3 - 0.035a; 4 - 4.22a; 6 - 
0.0001a; 7 + 19.08a; 8 + 0.21a; 10 + 148.844a; 12 + ...). 




(18) 



T = tOJ 



(19) 
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c). parallel /5-sheet ((1/A) sp =85.49 cm" 1 ) : p = 0.0028; = 0) < (for 

a; < 0.3 the expansion is Y ^ l s ^f l ~ 0.5a; 2 — 0.3a; 3 + 2.3a; 4 + 3.x 5 — 15.9a; 6 — 
14.6a; 7 + 57.65x 8 + 52.4a; 9 + CL64a; 10 - 158.2a; 11 + 403.3a; 12 + ...). 

These results mean the following. For all three structures the coupling interac- 
tion is surprisingly weak providing good conditions for the existence of a DB. 
For the a-helix and the anti-parallel /3-sheet the effective on-site potential is 
hard (goes more steep than a harmonic one) and the coupling interaction is 
repulsive (a non-zero value of a peptide group displacement tends to increase 
the values of the neighboring peptide groups displacements with the opposite 
sign). In accordance with the theorems of the paper [26] a stable DB must be 
of zig-zag shape (peptide groups undergo out of phase oscillations). For the 
parallel /9-sheet the effective on-site potential is soft (goes less steep than a 
harmonic one) and the coupling interaction is attractive (a non-zero value of 
a peptide group displacement tends to increase the values of the neighboring 
peptide groups displacements with the same sign). In accordance with the 
theorems of the paper [26] a stable DB must also be of zig-zag shape. Thus 
for no type of protein secondary structure a bell shape DB (peptide groups 
undergo in phase oscillations) can be stable. 

In what follows we restrict ourselves by considering the case of a DB in the 
a-helix because only for it the experimental data suggest such a phenomenon 
[11] (see the beginning of the Discussion). The set of difference equations (5) 
can not be solved analytically and is analyzed by direct numerical integration 
for 25 coupled equations (motivated by the fact that 25 is the average number 
of peptide groups in a-helicies of bacteriorhodopsin used in the experiment 
[11]). We use fixed-end boundary conditions X\(t) = a; 25 (t) = that seem to 
be more appropriate for the situation in proteins than commonly used periodic 
ones. We find as an example a 3-site DB. To prove that a localized and stable 
object can be obtained in the a-helix we resort to the usual procedure [25]. We 
define the discrete per site energy in 3 lattice sites around the 13-th central 
one 

1 14 

E(3) = ~ £ ei (21) 

1=12 

where e\ is the individual symmetrized local site energy density 

ei = \{^f + + (tf Xl ) + U(x t] x l+l )) (22) 

The existence of a localized and stable zig-zag shape DB is exemplified in 
Fig. 4 by the fact that the energy E(3) remains essentially constant with time 
and in Fig. 5 by an explicit picture. It should be stressed that for our realistic 
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potentials it is rather difficult to switch out the coupling interaction without 
modifying the effective on-site potential. Thus the commonly accepted way to 
construct a DB starting from an uncoupled (or else anticontinuous) limit is 
inapplicable in our case. That is why, instead, we find the initial configuration 
to fall into a DB from the very beginning at a given coupling interaction for 
the a-helix calculated in Sec. 3. In Fig. 4 and Fig.5 this configuration is chosen 
to be x 13 (0) = 0.3; xi 2 (0) = xi 4 (0) = -0.025 while ^(0) = for all other i 
and zero velocity at each site. The chosen amplitude x 13 = 0.3 corresponds to 
the deviations of the torsional angles A(p 13 = Aip 12 = x 13 /2 ps 10° from their 
equilibrium values. 



6 Discussion 

The results obtained testify that a localized and stable zig-zag shape DB can 
exist in all three types of protein secondary structure. However experimental 
data which can be interpreted with the help of the DB concept are available 
only for a-helix proteins but not /9-sheet ones [10], [11]. In our opinion the 
reason for this is in the way of excitation of long-lived oscillations in this 
experiments rather than in the propensities of different types of protein sec- 
ondary structure to sustaining DBs. The absorption of far IR laser pulse used 
in [11] by proteins requires long one-dimensional hydrogen bonded chains [17] 
that take place in a-helicies but are absent in /5-sheets. 

The revealed hardness of the effective on-site potential for the a-helix and the 
anti-parallel /9-sheet is in contrast with a widely spread point of view that 
such potential in biomolecules is mainly determined by hydrogen bonds and 
that is why should be a soft-type one [3], [26]. In our case detailed account of 
all interactions leads to a contrary result. 

The energy of the 3-site DB obtained for the a-helix is 3E(3)Iuj 2 ps 2 • 
10~ 12 erg m 40/csT at T=300 K° . This value means that at room temper- 
atures large amplitude DBs (Aip = Aip ps 10°) can hardly be spontaneously 
excited by thermal fluctuations. However the latter may not be true for DBs 
with smaller amplitudes and here it seems timely to touch upon the question 
of influence of temperature on the existence of DBs in protein secondary struc- 
tures. DBs are known to remain essentially stationary and localized and to be 
very long-lived at zero temperature. At non-zero temperatures thermal fluc- 
tuations tend to lead to DB motion and more rapid decay [2], [3], [33]. Within 
the framework of our interpretation of the the experiment of [11] with the help 
of a DB concept it should be recalled that the observed there long-lived os- 
cillations in a-helicies of bacteriorhodopsin excited by far IR laser pulse have 
a life-time of order of 500 ps accommodating approximately 1500 vibrations. 
The decay of a DB can be formally modeled by introducing a dissipative term 
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~^ dX df 1 * n ^ e hand side of the equation of motion (5) and corresponding 
white noise Q(t) in the right hand one. Here the friction coefficient 7 must be 
related with the white noise correlation function via the fluctuation-dissipation 
theorem < C»(*)Ci(0) >= 2kBTj6(t)5ij where ks is the Boltzman constant, T 
is the temperature, S(t) is Dirac 5-function and Sij is the Kronecker symbol. 
However such formal way of introducing dissipation does not reveal the micro- 
scopic origin of the friction coefficient 7. To do it for the case of peptide chain 
dynamics is a difficult task which is out of the scope of this chapter. Here 
we can afford ourselves only some preliminary speculations on this subject. 
There are different contributions into the mechanism of friction that can be 
conceivable for the peptide chain rotational dynamics in protein interior. The 
most obvious hydrodynamic damping by interaction with solvent seems to be 
of minor importance because at considered deviations of the torsional angles 
< 10° the linear displacements are much less than the diameter of a water 
molecule (see also [11] and refs. therein). In our opinion the dominant contri- 
bution into dissipation originates from the fact that the rotation of the peptide 
group at angle x round the axis a requires some translational displacement 
of axis points where the atoms iVj and C i _ l are initially situated (see Fig. 2) 
that leads to transvers distortions of the backbone (actually this is the excita- 
tion of transvers linear phonon modes). Such process consumes energy leading 
to an acoustic damping mechanism. At a macroscopic level this dissipation 
mechanism can be described by the so called "Landau damping" 7 ac oc uj 2 
[34]. However such phenomenological description requires support from mi- 
croscopic consideration that is beyond our knowledge at present. Moreover 
the situation with dissipation is not so unambiguous as one can conceive from 
above. There are testimonies that in some cases of nonlinear lattices the en- 
ergy relaxation from a DB first accelerates with the increase of the friction 
coefficient from zero value but then becomes slower and slower [5]! It is not 
clear at present whether this surprising result is generic or it is model depen- 
dent (e.g., dependent on the type of the thermostat used). In any case it opens 
interesting perspectives for application of the DB concept to the problem of 
localized vibrational energy storage in enzymes. In this connection the coinci- 
dence of the obtained value for the DB energy (rs 40/cbT at room temperature) 
with usual activation energies of enzymatic reactions is rather intriguing. One 
can imagine oneself that in secondary structures of enzymes there DBs can 
be excited at substrate binding and be long-lived enough to play some func- 
tional role as the sources of energy localization and storage utilizing in such 
a way the binding energy. For instance a DB may serve as the so called "rate 
promoting vibration" (RPV) for enzymatic reactions involving enviromentally 
assisted hydrogen tunneling. The RPV acquires more and more experimental 
testimonies [35], [36], [37], [38], [39] but its physical origin still remains mys- 
terious. However the estimated dominant peaks in the spectral densities of 
the RPV indicate motions on the 150 cm~ l frequency scale [38] that is rather 
close to the obtained frequency 115 cm~ x of the DB in the a-helix. Here we 
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refrain from further discussing this speculation. 



7 Conclusion 

The conclusions of this chapter are summarized as follows. The dynamics of a 
peptide chain is considered within a realistic model with stringent microscop- 
ically derived coupling interaction potential and effective on-site potential. 
The coupling interaction is found to be surprisingly weak for all three main 
types of protein secondary structure but different in character: repulsive for 
a-helix and anti-parallel /3-sheet structures and attractive for parallel /3-sheet 
structure. The effective on-site potential is found to be a hard one for a-helix 
and anti-parallel /3-sheet and a soft one for parallel /3-sheet. In all three types 
of protein secondary structure a stable zig-zag shape discrete breather asso- 
ciated with the oscillations of torsional (dihedral) angles can exist due to the 
weakness of the coupling interaction. 
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Fig. 1. A schematic representation of the peptide chain with all designations neces- 
sary for the text. The covalent bonds are shown as bold lines. Dotted lines denote 
hydrogen bonds. Dashed lines denote auxiliary axes. 
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Fig. 3. The normalized effective on-site potential V e ff(x) (see (3)) for the a-helix. 
Here I is the moment of inertia of the peptide group relative the axis a and to is the 
frequency. The range for the angular displacement x (defined in (1)) is 4ir because 
it is twice of the torsional angles tp and tp (see (1)). The part of the potential over 
the upper cut off has a maximum with the value « 1.3. 
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Fig. 4. The dependence of the energy in the 3-site discrete breather determined by 
(21) and (22) on dimensionless time r (19) in the log-time scale for the a-helix. 
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Fig. 5. The explicit picture of the long time (in dimensionless units (19)) behavior 
of the discrete breather in the a-helix of 25 peptide groups named here as sites. The 
displacement x (defined in (1)) is the angular deviation of the peptide group from 
its equilibrium position. 
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